NONLINEAR ELASTICITY FROM ATOMISTIC MECHANICS 



C. ORTNER AND F. THEIL 



Abstract. We present sharp convergence results for the Cauchy-Born approximation of gen- 
eral classical atomistic interactions, for both static and dynamic problems, for small data. 



1. Introduction 

The Cauchy-Born model (or, approximation) for Bravais lattices is the most widely used 
nonlinear elasticity model of crystal elasticity. It is obtained from the Cauchy-Born rule: the 
stored energy per unit volume under a macroscopically homogeneous deformation equals the 
energy per unit volume in the corresponding homogeneous crystal. In some simple cases the 
Cauchy-Born approximation for Bravais lattices can be justified as a statement about minimum 
energy states [HIE]. This means that microstructural relaxation effects are ruled out. In less 
restricted settings it has been shown to provide a highly accurate approximation of crystal 
elastostatics [TU] . 

It is highly desirable to obtain comparable results for evolutionary problems because of the 
importance of phenomena linked to energy transport and dissipation in crystallographic lattices 
[1'6\ [TBI [12] • The main contribution of this article is a rigorous approximation error analysis 
of the Cauchy-Born wave equation compared to the Newtonian equations of motion in the 
atomistic model. In this pursuit we are inspired by a recent effort of Blanc, Le Bris and Lions 
[3] who propose this problem and solve it for various one-dimensional examples with different 
classes of pair interactions. 

Our own convergence result is formulated for a general class of multi-body interactions in 
an infinite lattice, and only requires the assumptions that the "reference lattice" is a "stable" 
Bravais lattice (cf. Section 15. 2p and that the interaction strength decays sufficiently fast. We 
also provide a rigorous approximation error analysis of the Cauchy-Born approximation for 
static problems and the same conditions. The lattice stability assumption is essential and 
cannot be removed (see § [C|) . 

Our main dynamic result. Theorem [621 is concerned with the Newtonian equations of motion 
and the nonlinear Cauchy-Born wave equation. In the scaling limit where both spatial and 
temporal fluctuation are of the order e (the lattice spacing) we prove that the atomistic solution 
converges a solution of the Cauchy-Born wave equation. Moreover, difference between the 
atomistic and the Cauchy-Born solution is of order e^. The result for the static case (Theorem 
15. 4p is analogous. 

These strong results are the consequence of sharp quantitative links between particle models 
and continuum models, which are also useful in a wider context |28[ \20\ I29j . In particular 
the localization technique ()1.2p provides a continuous representation of discrete objects. A 
key concept is the notion of atomistic stress, which gives a natural weak form of the first 
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variation. Pointwise second-order consistency of the Piola-Kirchhoff stress of the Cauchy- 
Born model with the atomistic stress is established in Theorem 14.31 Error estimates between 
local minimizers of the Cauchy-Born and atomistic model are an immediate consequence of 
consistency. The proof in the dynamic case is based on a similar pointwise consistency result 
for the divergence of the stresses. 

1.1. Main results. We assume that the scaled total atomistic energy can be written as 

where A := Z'^, n : eA — )■ M"^ is a discrete displacement and V a multi-body potential describing 
the interaction between a site and the rest of the body. The continuum limit is characterized by 
the Cauchy-Born energy density function W{f) = ^({F/?}p£A\{o}) the associated functional 
E{u) = Jjg<j W^Vu) dx. We will show that E characterizes the asymptotic behavior of solutions 
of static and dynamic problems associated with E'^ . We assume throughout this section that 
d < 3 and that A is stable (see 



Theorem A (Elastostatics). There exists constants Csta.t,C,So > such that for all f : 
^R"^, f : eA^ satisfying 



\i<So and <Cstate , 

and for e sufficiently small, there exist local minimizers u, u'^ of the energies 

E{u)- f f-udx and E' {u) - e'' f {£) ■ {i) , 



which also satisfy the bound ||Vn — Vn'^||^2 < Ce^ . 

The definition of the norms || • ||i, || • ||-i and a precise formulation of the result are given in 
Theorem [Ol 

Theorem B (Elastodynamics). There exist constants n,T,C > such that, for all initial 
states {uq,ui), (nQ,n|) satisfying the bounds 

\\^uo\\l°° < K and ||Vuo - Vno||^2 -h ||ni - 11^2 < Qnit e^, 

where Cinit is independent of e, and for all sufficiently small e, there exist unique solutions 
u € C2([0,T],if2(R'i)) and G C^{[0,T], f{A)) of the Cauchy problems 

a - divS'=(Vn) = 0, u{t = 0) = uq, u{t = 0) = ui, 

ii^ + SE'{u^) = 0, u%t = 0) = u'q, u%t = 0) = uf , 

where 6E^ denotes first variation of E^, and S^j{F) = ^^(F) is the first Piola-Kirchhoff stress 
tensor in the Cauchy-Born model. The solutions n, satisfy the estimate 

^max (|[Vn(t) - Vu^(i)||i2 + \\u - ^^(OIIlO < Ce^. (1.1) 

The definition of the norm || • ||3 and the precise formulation of the result is given in Theo- 
rem [621 We note that the condition ||Vno||Loo < k is a fairly mild condition, which only ensure 
that uq is a "stable configuration". 

We make no statement about the maximal time interval for which the estimate (jl.ip holds. 
Such a statement could be made provided one establishes an atomistic Carding inequality. 
Moreover, such a result would also allow us to treat large deformations in the static case. We 
stress, however, that our estimates hold for a macroscopic time interval. 
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1.2. New Concepts. The comparison of discrete displacements with continuous displace- 
ments u is achieved by the introduction of several approximation operators derived from a 
nodal basis function ( with compact support (see ^ l2.1l for the details). We define the Lipschitz 
continuous interpolation and a quasi- interpolation u^, which has a Lipschitz-continuous gra- 
dient. The key property property of the framework is that it delivers an integral representation 
of finite differences 

v{^ + ep)-v{0 ^ I ^1 (^)v,t;(x) dx, where x\ := e"^ /' C(e + tep - x) dt. 

£ JW^ ' Jo 

(1.2) 

The kernel x| p can be understood as a mollified version of the line measure supported on the 
bond {C + tp\t e [0,e]}. 

A version of this technique was proposed by Shapeev [31] for the construction of atomistic- 
to-continuum coupling schemes. Here we propose an analytically convenient and stable variant 
of his idea, which is analyzed in detail in [25]. 

Using ()1.2p we can construct an integral represenation of the first variation of the atomistic 
energy, 

{5E'{u''),v) = [ S"(n=;x) : V?;(x)dx, (1.3) 

for all compactly supported virtual displacements v, where is the stress associated with a 
discrete displacement u^, 

S%u';x) = E E ® - OUlpix); (1.4) 

?eeAeeeA\{5} 

here, /^^g' is the force acting between atoms ^ and ^' due to the site-energy associated with 
atom ^. A precise definition of the atomistic stress is given in § 14.11 We will prove in § 14.21 
that IS'^ - S^l = 0{e^) and in § [63] that |div(S= - S=)| = ©(e^), which are the key technical 
ingredient required for proving Theorems A and B. 

Remark 1.1. The tensor is closely related to the stress defined by Hardy [12] . Indeed, 
equation (4.5) in that work is essentially an Eulerian version of ()1.4p . for pair interactions, with 
a generic weighting function x| p- A general account of stress in molecular dynamics simulations 
is given in the recent work of Admal and Tadmor [1]. 

A concern discussed in [1] is the non-uniqueness of stress: note that is defined by ()1.3p 
only up to a divergence-free tensor. Indeed, it is possible to add "discrete null-Lagrangians" to 
the atomistic energy, or decompose the atomistic energy in different ways, which would lead to 
different definitions of the atomistic stress . 

Concerning this question, our work provides a selection mechanism based on comparison 
with the Cauchy-Born stress. If V is selected to be "as local as possible" (cf. § I2.3.4p and 
to satisfy an inversion symmetry (cf. § I2.3.3p . then the atomistic stress S*" and the Cauchy- 
Born stress S'^ are second-order close. Since there is a natural definition for S'^, our version of 
the atomistic stress is reasonable whenever the atomistic configuration is "close" to an elastic 
continuum configuration. □ 

1.3. Outline. Since we work in an infinite domain, and admit an infinite interaction range, 
even the definition of the atomistic energy is non-trivial. Section [2] is devoted to this task. At 
the same time we establish various auxiliary results that are useful for the subsequent analysis. 
In Section [3] we define and analyze the Cauchy-Born approximation; in particular, we establish 
differentiability of the stored energy function and establish a convenient functional analytic 
setting. In Section [5] we derive and analyze the atomistic stress, which plays a prominent role 
in our analysis. In Section [5] we present a rigorous approximation error analysis of the static 
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Cauchy-Born approximation. Finally, in Section [6] we establish approximation error estimates 
between the solutions of Newton's equations of motion and the Cauchy-Born wave equation. 

1.4. Summary of notation. Throughout, M denotes the real numbers, and d € N the space 
dimension. 

Positions in space are usually denoted by x,y, z E W^, while lattice sites are denoted by 
^,7] G A := Z'^. Displacements, either continuous or discrete, are denoted by u,v,w. Lattice 
directions are denoted by p,<^,T G A* := Z'^ \ {0}. 

Matrices are denoted by capital letters, A, F, G, S € M'^^'^. In particular, A is reserved for the 
lattice orientation (beginning of §[2]), and S for stress tensors. 

If / : M'' ^ is differentiable in X then we denote its Jacobi matrix by Vf{x) and a 
direction derivative by Vrf{x) = Vf{x) -r, r £ W^. Higher derivatives are denoted by f and 
are understood as j-linear forms with range in M™. 

If £/, M are Banach spaces and ^ : — )• ^ is j times Frechet differentiable (or sim- 
ply, differentiable) in a then we denote its j-th derivative (or, variation if ^ = M) by ^ ^ 
which is understood as a j-linear form with range in SS. If a, oi, . . . , Oj G =(2/, then we write 

^ {a)\a\^ . . . , Oj] to evaluate this form. 

If £ : ^ ^ R is a linear functional, then we write 1(a) = (£, a). If E : £/ — > M is differentiable, 
then we write SE{a)[ai] = {5E{a),ai). If it is twice differentiable then we will also write 
6'^E{a)[ai,a2] = {6'^E{a)ai,a2). 

We use LP, W^'^ , for p G [1, oo], to denote the standard Sobolev and Lebesgue spaces, usually 
on the domain R'^ (and otherwise specified). We will also employ the so-called homogeneous 
Sobolev spaces, W^'^, which are defined in § 13.11 Negative-norm (dual) spaces are denoted by 
W-'^'P = {W'^'P')* and W-'^'P = {W^'P')*, where p' := p/{p - 1) G [l,oo] is the dual Sobolev 
index. 

The discrete "Lebesgue space" is denoted by ^P , usually with domain and otherwise 
specified. Since most discrete versions of W^'P would be equivalent to iP, we will not define 
such a space. However, we make heavy use of discrete homogeneous Sobolev spaces W^'P, which 
are defined in ^2.11 

Finally, we make the convention that "<" stands for "< C" , where C is a generic constant 
that may not depend on any data in the model, nor on any functions involved in the inequality. 
In particular, we will make explicit all dependence on the interaction potential, which is crucial 
since we admit an infinite interaction range. 

2. The Atomistic Energy 

We formulate an atomistic model with classical multi-body interactions on an infinite Bravais 
lattice. As discrete reference domain, we choose A := l/, where (i G N (later restricted to d < 3) 
is fixed throughout. We admit deformations of the form y{^) = fK^+u{^), where u is an unknown 
displacement and A G R'^^'^, det A > 0, defines the reference state of the system, A • A, which 
may be an arbitrary Bravais lattice. 

We now present a formal definition of the atomistic potential energy, which we make rigorous 
throughout the remainder of the section. Let A* := A \ {0} denote the set of lattice directions. 
For discrete maps f : A — )■ R"*, m G N, we define the finite differences and finite difference 
stencils 

Dpv{i) ■.= v{^ + p) -v{i), for^GA,pGA*, and 
Dv{0--={Dpv{0)p^K,, for^GA. 
A convenient space of finite-difference stencils is 

& -.= {9 = {gp)peA, I gp e R"^, \\gy < oo}. 
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equipped with the norm \\g\\c^ := maXpgA* |ffp|/|/o|- 

Next, we assume that there exists a site-energy F : ^ — t- M U {±00}, so that the atomistic 
potential energy of a displacement n : A — > can be written, formally, as 

E''{u) := ^<^^{u), where ^^(n) := V{Du{^)). (2.1) 

Note that the site-energies are not well-defined for general u, and moreover their sum need not 
exist. In the remainder of this section, we introduce a discrete function space setting in which 
we can make (|2.ip rigorous. 

We remark that V implicitly depends on A, but since A is fixed throughout, we suppress this 
dependence. 

2.1. Interpolation of lattice functions. We denote the set of all vector-valued lattice func- 
tions by W and those with compact support by Wq: 

W := {v : W^} and := {v e W \ supp(i;) is compact}. 

To facilitate the transition between continuous and discrete maps we introduce two (quasi-) 
interpolants of lattice functions. 

Let C G W^'°°{R'^; M) satisfy C(6 = for ^ G A \ {0} and C(0) = 1. We understand C( • - 
as a nodal basis function associated with the site ^, and define the first- order interpolant 

v{x) ■.= ^v{^)C{x -C), iorveW. (2.2) 

?6A 

We shall assume throughout that C ^ 0? C has compact support, C is symmetric about the 
origin Ci~x) = C{^)i a^id that the associated interpolation operator reproduces affine functions: 
^^g^(a-|-6-^)C(x — ^) = a + b-x for all a G M, 6 G W^. The latter property implies, in particular, 
that Jjgtj C dx = 1. 

Following [31j . we also define a quasi-interpolant obtained through convolution of v with 

v{x) := {C * v){x) = C{x-y)v{y)dy. (2.3) 

In general, w(^) 7^ v{^) for ^ G A, hence it is only a quasi-interpolant. The introduction of this 
second interpolant leads to a set of techniques centered around the localization formula \61\ [25] 
(cf. the scaled version in (|1.2|) ) 

^P^^(0 = / Xc,p(a^)Vpv(x)d2;, where x^A^) ■= [ + tp - x) dt, (2.4) 

JRrf Jo 
which yield surprisingly strong connections between the atomistic model and its continuum 
counterpart (in particular, the definition and analysis of the atomistic stress in To prove 
(|2.4p . we simply note that 

Dpv{i)= f Vpv{^ + tp)dt= [ [ C{^ + tp-x)dt\/pv{x)dx. 

Jo JR'^ Jo 

Remark 2.1. A canonical choice for is the Ql-nodal basis function 

Cix) :=ntimax(0,l-|x,|), 

then supp(C) = [—1, l]'^ and C is piecewise multi-linear. In this case {v\v £ W} is the Ql finite 
element space, or equivalently, the space of tensor product linear B-splines and {v\v (z is 
the space of cubic tensor product B-splines |14j . However, other choices are equally possible, 
and indeed necessary in some situations |28] . Since none of our results require explicit knowledge 
of the type of interpolant, we admit the most general case. □ 
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We now collect several auxiliary results on the lattice interpolants introduced above, all of 
which are established in |25| . 

Lemma 2.2. Let v € W, then its first-order interpolant ()2.2p belongs to W^^^{W^;M.'^), and 
V € 14/j^'~(M'^;M'^). Moreover, for any p € [1, oo], 

ll'^^ll^p ^ W^Wip ^ W^Wlp < II^^IIlp ^ Ibll^p ^"^^ € W , and (2.5) 
IIVuIIlp < ||Vn||LP < ||Vf IIlp Vt; G (2.6) 
(All constants in the above estimates are independent of p.) 

Next, we state a useful embedding result; of particular interest is the case j = m = 1 and 
q = cx), which states (employing (|2.6|) ) that ||Vn||ioo < ||Vti||LP for all p G [l,oo]. The proof 
uses the fact that the "mesh size" in Z*^ is one, and that £p C 

Lemma 2.3. Let < j < m < 3 and p < g G [l,oo], then ||V™m||l<j < ||V''u||iP for all 

2.2. The space of admissible displacements. Since the atomistic model is formally trans- 
lation invariant, we define equivalence classes 

[u] := {u + t\t eM'^}, foiueW, 

and, for p G [1, oo], define corresponding function spaces 

W^'P := {[u]\ue W, \\Vu\\lp < +oo}. (2.7) 

We will not make the distinction between u and [u], whenever it is possible to do so without 
confusion, for example, when a statement or function with argument u is translation invariant. 

Proposition 2.4. 'W^'P is a Banach space. For p G (l,oo), the subspace {[v] \ v G Wq} is 
dense in W^'^. 

Let uG#'^'^,j9<oo, bea discrete displacement and y(^) := + m(0 be the associated de- 
formation. It is shown in [25] that \u{(,)\ <C |^| as |^| — >■ oo. Therefore, the discrete deformation 
satisfies the far-field boundary condition 

y(e) = A^ + o(|e|) asl^l^oo. 

Aside from satisfying a boundary condition we also require that deformations are injective. 
Since our main results do not cover arbitrarily large deformations, we will circumvent the 
question of injectivity by placing an L°°-bound on the displacement gradient (this will be 
ensured through conditions on the external forces). Notationally, we fix a constant k > 0, and 
define 

:= {ueW \ \Dpu{S,)\ < K for ah ^ G A,p G A*}- (2.8) 
If K is chosen sufficiently small, then displacements belonging to give rise to injective 
deformations. Indeed, if k < ||A~-^||~-^, where || • || denotes the ^^-operator norm, then 

\yiO - yiv)\ > i^\H^ - v)\ v^,77gA, (2.9) 

where := 1 — k||A^^|| > 0. 

To conclude the discussion of discrete function spaces, we note that we can define a smooth 
nodal interpolant on z^^'^, which will be useful in interpreting our results. 

Lemma 2.5. Let u G W^'^, p G [l,oo], then there exists w G W^'^ such that u = w =: lu. 
Moreover, 

I|Vu||lp < \\VIu\\lp < IIVmIIlp Vn G W^'^. (2.10) 



NONLINEAR ELASTICITY FROM ATOMISTIC MECHANICS 



7 



2.3. Assumptions on the interaction potential. 

2.3.1. Energy difference. We assume throughout that V^(0) = 0. Physically, this condition 
means that is an energy difference between the deformed state = A,^ + u{S,) and the 
reference state y{£,) = A^, which may have infinite energy. 

2.3.2. Smoothness. We assume that V is "smooth" at injective configurations; more precisely, 
we define '•= {g ^ ^ \ 11911^ < and assume that V G C^{S!^), for some k>2 and for all 
K that are sufficiently small. 

For g € and for any "multi-index" p G Ai, 1 < j < A;, the partial derivative 



dgpi ■ ■ ■ dgp^ 

exists; Vp{g) is understood as a multilinear form acting on families of vectors h = {hi, . . . , hj); 
Vp{g)[h]=Vp,...p^{g)[hi,...,h,]. 

We also define the associated partial derivatives of the site-energies by ^^^p{u) := Vp{Du{^)), 



or, 



where -Dpf (0 = {Dp^vi{0, ■ ■ ■^Dp^Vj{i)) for pi e A*, G #, z = 1, . . . ,j. 

2.3.3. Symmetry. We assume throughout that V satisfies the following inversion symmetry: 

F((-<7-p)pgA.) = V{g) yg e ^. (2.11) 

This condition is physically motivated by the fact that permutations of atoms and isometrics 
should not change the energy of a system. The requirement ()2.1ip then assumes that the global 
energy was partitioned in a way that preserves this symmetry. 

Lemma 2.6. Let F G M'^^'^, |F| < k, then 

y_p(F- A,) = (-iyVp(F-A,) VpGAi, 1 < j < k. (2.12) 



Proof. Let g^ := F- A*, then g'p := (— F(— p))pgA, = fl'F- Since V is k times differentiable in 
we can differentiate (|2.1ip at gf. Evaluating the first and second derivatives gives the stated 
result. □ 

2.3.4. Decay Hypothesis. Sufficiently rapid decay of the interatomic interaction is a crucial 
ingredient in our analysis. We define the following basic bounds on the interaction potential: 
With this notation, we define the bounds 

"i(p) —DLi sup ||Vp(si)|| forpGA^, l<j<k, 

where := sup^g^jgd^j ^[/ii, . . . , /ij] for a j-linear form i. We shall assume 

throughout that V and k are such, that 

^(j) ^(p) < foj. i<j<k, (2.13) 

which will ensure that E'^ is k times Prechet differentiable (cf . Theorem 12. 8p . 

In order to describe the class of admissible potentials, we also discuss the decay assumption 
we will require in the static and dynamic approximation error analysis. To simplify the following 
expressions, we define \p\ := "^l^i \pi\, for p € Ai, j € N. 
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Let p G [1, oo] and 2 < j < k. In the static analysis, we will assume finiteness of certain 

M^^'P^ := rniP\p), where 
pGAi 

1 

m(^)(p) ■.= mip)\p\'^{j:l=2 (IPi X P^\ + \Pi\ + for p G A^. 

In the dynamic analysis we will assume finiteness of certain constants of the form 
mU,p) ._ rn^^\p), where 
peAi 

^f^P) ■■='^^^^{n=2 {\Pi X P^\ + IPil + Ip.I))^, for P e Ai 

These constants naturally arise in the modeling error estimates established in ^4.31 and 

We stress that, while finiteness of M^^\ 1 < J < ^, is a standing assumption, we will assume 
finiteness of M^^'^^ and M^J'^\ for certain choices of j and p, only when required. 

Remark 2.7. If contains only pair interactions, then we can write V in the form 

peA, 

which clearly satisfies the symmetry (|2.1ip . Moreover, we show in ^B.ll that. if ip^^\r) < r~°'~^ 
for r > 1, 1 < J < A;, then 

pGA. pgA* 

It therefore follows that the constants M^^'^ are finite provided that a > d and that Mi"''^^ M^'^^ 
are finite provided that a > d+5/2. In particular, this implies that the Lennard-Jones potential, 
</?(r) = r~^^ — 2r~^, is included in our analysis. 

In ^B] we discuss other commonly employed potentials and show that they fall within our 
assumptions. □ 

2.4. Definition of tlie atomistic energy. We mentioned at the beginning of Sj2]that, due 
to the infinite domain and the infinite interaction range, the definition of the energy (j2.ip is 
non-trivial. The purpose of this section is to give ()2.ip a rigorous interpretation. 



Theorem 2.8 (Regularity of E^). (i) If u e Wq n , then ($g(n))5gA G ; that is, E^( 



u 



given by (|2.ip is well-defined. 

(a) E^ : {Wq n J^, II • 115^1,2) — 7> M is continuous; that is, there exists a unique continuous 
extension to W^"^ n .J^ , which we still denote by E^. 

(Hi) G C'=(#i'2 n JT), with 

b^E'^iuM = E E ^U^)[DpVm (2.14) 

?gA pgAi 

for j = 1, . . . , A:, for all u G W^'"^ and v G #0^. Moreover, if 1 < j < k and 'Ylil=i ^ — 
Pi G [1, 00], then 

j 

\6^E^{u)[v]\ < M^^^HWVviWLP.. (2.15) 



i=l 
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The proof of Theorem 12.81 is estabhshed throughout the remainder of the section. We begin by 
estabhshing a simple bound on finite differences, which gives a first ghmpse of the locahsation 
technique used at crucial steps throughout the paper. 

Lemma 2.9. Let p € A* and p € [l,oo], then 

\\Dpu\\gp < \\Vpu\\LP < \p\\\\/u\\lp Vn e 



«6A 5 

< 



Proof. The case p = oo is trivial, hence suppose that p < oo. Using the localisation formula 
(|2.4p . the fact that {({ ■ — ^) | ^ S A} is a partition of unity (and hence — dx = 1), 

we can estimate 

JZl^p'^iOf = Y.\ [' I C(^ + tp-^)Vpn(x)dxd/ 
y] / / C,{i + tp- x)\V pu{x)\'' dxdt 

^ Jt=0 JR'i 

= [ \Vpu{x)\^dx <\p\P\\Vu\\l^. □ 
Lemma 2.10. Let A, and Vi G #1'^% 1 <i<j<k, with J2i=i ^ = 1' ^^e'^ 

^eApgAi 

Proof. Let "Wj G yf/^^Vi such that iDj = /fj, then 



i=l 



i=l 



Summing over ^ G A, p G Ai, applying the generalized Holder inequality, followed by Lemma [2T9l 
yields 



Applying ()2.10p yields the stated result. □ 

Remark 2.11. Lemma [2 . 1 U I implies . in particular, that we can interchange the summation in 
series of the form X^^gA X^pgA^ '^p{9£){^p'"{€)\^ provided that g^ G and Vvi have sufficient 
integrability. We will henceforth perform interchanges of summations without further comment. 

□ 

Lemma 2.12. Ifu^W^, then 

E • ^P<^) = 0- (2.16) 

geApeA, 
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Proof. According to Lemma l2.10l the sum on the right-hand side of (j2.16p is weh-defined, the 
summand belonging to i^. 

Since <^>^^p(0) = Vp{0), and interchanging the order of summation, we obtain 

^a^) ■ D,u{i) = Y V,{0)YDpu{0. 

CeApGA, pGA* 5gA 

For u € Wq the right-hand side clearly vanishes. □ 

Motivated by the previous lemma, we define 

i-(n) := Y ^?(«)' where ^^(n) := <^^{u) - ^ ^CpW " Dpu{^), (2.17) 
CgA peA, 

then we immediately obtain the following result. 

Lemma 2.13. € C''{W^'P n J^T), for all p G [1, oo], with 

^i,p{u) = %p(n) - $5,p(0), for p e A*, and 
^^,p{u) = <^^,p{u), for peAi,2<j<k. 

Proof. This result follows immediately from the fact that $g S C''(y^^'^n^) and the estimate 
Y \'^U0)\ \Dpu{0\ < Y lpk(p)l|V^x||LP = M«||Vn||iP. □ 

pGA, pGA, 

Lemma 2.14. Let u G n then ^^{u) € that is, (j2.17p is well-defined for all 

u G #^'2 n jr. Moreo?;er, G C^(#^'2 n mi/i variations given by ([231]) ("wii/i ^J'" 
replaced by E^). 

Proof. Fix u G #^'2 n^, then 6'Z)n(C) G for ah ^ G A and 6 G [0, 1]. Hence, we can expand 

|.g(n) = 1.5(0) + Y ^5,p(0) • DpuiO + Y Vp,{9c)[Dpu{0,DM0i (2-18) 

peA, p,fGA, 

where = 9^Du{i) G for some % G (0, 1). By definition, = <&5,p(0) = 0, and hence, 

p,?gA* 

It now follows from Lemma 12.91 that ^^(u) G that is, E^{u) is well-defined. 

The differentiability of $ can be shown using analogous arguments. □ 

Proof of Theorem [Ol (i) If n G ^ fl ^ then Lemma [2Tll implies that <I>g(n) G while 
Lemma [27T2] ensures that {^^{u) - ^^^{u)) G 

(ii) For u G PI J(f , Lemma 12.121 implies that E^{u) = E^{u). Since E^ is continuous on 
y^^'^ n J(f , and Wq is dense in W^'"^ , it follows that E'^ is the continuous extension of E^ to 

(iii) Since E'"^ = E^, this statement follows from Lemma 12.141 The simplified representation 
of 6E^ (replacing p with $g,p) follows from Lemma 12.121 For j > 1 and p G Ai we have 
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Remark 2.15. 1. From the proof of Lemma 12.141 is becomes apparent why we assumed in 
^2.3.21 that k > 2. Only in that case are able to show that is well-defined on W^'"^ fl J^. 
We note, however, that if we only assume k = 1, then the the definition (j2.ip still yields 

2. To define E^{u) it is not necessary to assume a bound on ||Vn||/,oo. Indeed, if we require 
that V is smooth at all configurations Du{S,) for displacements u satisfying the injectivity 
requirement ()2.9p . then Theorem 12.81 can be extended to this class of deformations. 

The key observation is that |p|~"^|L'pti(^)| — )> uniformly in p as |^| oo, and also 
|p|~^|I5pu(.^)| uniformly in ^ as \p\ — ?> oo. (The first statement follows from the fact 
that ^ 1-^ Dpu{^) G the second statement follows from the inequality |/9|~^|-Dpu(^)| < 
|/9|~^/^||Vu||£,2, which is easily established from ()2.4p .) In particular, this implies that Du{^) € 
Si^ for |,^| sufficiently large. One can now apply the expansion ()2.18p in the far-field. □ 

3. The Cauchy-Born Approximation 

The Cauchy-Born elastic energy density function W : W^^'^ — ?• M U {±00} is defined by 

W{?) K). (3.1) 

In the regime of "smooth elastic" deformations, the Cauchy-Born model is a popular approxi- 
mation to the atomistic model E'^ [3l fTOl [321 fT5] . 

Formally, if u : M"' ^ R'^ is smooth, then V{Du{C)) ^ V{Vu{C) ■ K) = W{Vu{C)) and hence 
we can approximate E^{u) by 

E%u) := [ W{Vu)dx. (3.2) 

In the remainder of this section we introduce a function space setting to make ()3.2p rigorous, 
and establish associated auxiliary results. 

3.1. Homogeneous Sobolev spaces. The formal Euler-Lagrange equation associated with 
E'^ is a second-order elliptic system in M'^. Due to the translation invariance, convenient function 
spaces for equations of this type are the homogeneous Sobolev spaces (or, Beppo-Levi spaces [8]). 
Here, we use a variant of spaces of equivalence classes: 

W""'^ ■={[u]\ue W™'^(M'^;M'^),Vu G W"^"^'^}, for m = 1, 2, . . . . 
The space VF*"'^ is equipped with the norm 

It is easy to see (see [27] for a proof of the case m = 1; the general case is analogous) that 
W"^'^ is a Banach space, and that the subspace {[w] |f G C^} is dense, where := {u G 
C°°(]R'^;R'^) |supp(n) compact}. 

The natural space of continuous displacements is W^'"^. In order to avoid non-interpenetration 
of matter we shall assume that all displacement gradients satisfy a uniform bound. To that 
end we define 

K ■.= {ue ||Vn||L°o < k}, 

where k is the same constant as in the definition of (|2.8p . 

A straightforward extension of \27\ Thm. 2.2] shows that |u(rE)| <^ \x\ as |x| — >■ 00 for all 
u G W^'^ n K, hence the associated deformation y{x) = Ax + u{x) again satisfies the far-field 
boundary condition y{x) ~ Ax as \x\ — )■ 00. 
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3.2. Definition of the Cauchy— Born energy. In this section we make the definition of 
the Cauchy-Born energy ()3.2p rigorous. We first analyze the stored energy function W of the 
Cauchy-Born model. 

Lemma 3.1. W e C''{{f e M'^^'^||F| < k}). For |F| < k and (Gi,...,Gj) G {R'^'^^y , 
1 ^ J < we have 

6miF)[Gi,...,G,]= ^p(F-A,)[Gipi,...,G,pj]. 

peA-i 

Proof. This result follows immediately from the fact that, for |F| < k, F • A* € D 
The derivative of W can be represented by the first Piola-Kirchhoff stress tensor, 

for i, Q = 1 . . . d, u & K. (3.3) 

F=Wu{x) 

Then we can write (5Vl^(Vii)[G] = S'^(n) : G, where ':' denotes the usual Frobenius inner product. 

Following the arguments in H2A\ we obtain that is well-defined in H K and has a 
continuous extension to W^''^ fl K. The proof is analogous to the proof of Theorem 12.81 

Proposition 3.2 (Definition of E"). (i) IfueC^nK, then W{Vu) G L^; hence E''{u) is 
well-defined by (j3.2p . 

(ii) E'^ : (Cq° n iT, II • ||^yi,2) R is continuous; hence there exists a unique continuous 
extension to W^'^ H K, which we still denote by E^. 
(Hi) E"^ eCHW^'^nK;\\ ■ + ' W^i^oo) with 

S^E^iuyv] = [ 6m{\/u)[\7vi, Vvj] dx, for v G {C^y , I < j < k. (3.4) 
Moreover, for "^j—i ^ = 1, we have 

j 

6^E%u)[v] < m(^') Yl llVuillLPi. (3.5) 

i=l 



Sf„(u;x) := 



dW{F) 



4. Stress 

For u G W^'"^ n K the canonical representation of dE'^{u) is (cf. (|3.4|) with j = 1) 

{6E%u),v) = [ S"(u) : Vvdx, for veC^. (4.1) 

The first Piola-Kirchhoff stress is dual to the virtual displacement gradient Vu . By contrast, 
the first variation of the atomistic energy is more commonly expressed in terms of forces, which 
are dual to v. The purpose of this section is to derive and analyze an atomistic concept of 
stress that yields a representation of 5E^ analogous to (j4.ip . This will then be employed in a 
sharp consistency analysis of the Cauchy-Born approximation in § 14.2^ § 14. 3^ and § 16. 3i 
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4.1. An atomistic stress function. The "canonical weak form" of given in ()2.14p . is 

{5E^{u),v) =YY1 ■ Dp<0, for V G n- (4.2) 

5eA pGA, 

To proceed, we fix u G #5 but test 5E^{u) with v instead of v. We apply the localization 
formula (12.41) to obtain 



= / I E E [^CpW ^Hx5,p(^)| : Vt;dx. (4.3) 

Using the decay assumption (|2.13|) one can apply Fubini's theorem to justify the interchange 
of integral and sums, and thus obtain the following result. 

Proposition 4.1. Let u e n JT, then 

{6E^{u),v)=[ S^('u;x) : Vt;(x)dx Vu e ^o, 



where S''{u;x) -=^^2 [^«'p(^) P]xiA^)^ 

geApeA, 



(4.4) 



and x^,p is defined in ([2 

Remark 4-^- The fact that the test function in the left-hand side and right-hand side 
of ()4.4p differ may seem counter-intuitive at first. Allowing this seeming discrepancy makes the 
rather simple definition of atomistic stress possible, and will lead to sharp consistency estimates 
requiring only modest analytical effort while the subsequent error analysis in f|5]and ^requires 
only minor adjustments. □ 

4.2. Second-order accuracy of the Cauchy Born stress. In this section we prove the 
following pointwise second-order consistency estimate between the Cauchy-Born and atomistic 
stress functions. 



Theorem 4.3. Let u G W^^.oo p ^ ^ 



en 



|S>;x)-S^(n;x)| < ^ mi^\p) WV'yU^^,^,^ 



(4.5) 



P6A2 



where mi°°'^ is defined in ^2.!^.4\ and for each p, G A* the set Vp^ satisfies Up^ C -6(0, c(|/9|-|-|?|)), 
for some constant c > 0, — fp<^ = fp,;, and 



vol(i^pJ < |p X ?| IpI + |?| 



(4.6) 
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Before we embark on the proof of Theorem 14.31 we estabhsh two useful identities for the 
weights x^.p) which will enable us to exploit the inversion symmetry (j2.1ip . 

Lemma 4.4. Let x,/9 G ; then 

Ex?.p(^)(^-^)= -^P- (4.8) 
?6A 

Proof. Both results rely on the assumption that afhne functions are invariant under the first- 
order interpolant. Clearly, this is still true on a shifted grid: if affine, then for 
any z, a; € M'^ we have 

v{x) = C{x — r])v{r]). (4.9) 

r)G(A+2) 

To prove ()4.7p we write out the left-hand side, and employ (|4.9p with v{x) = 1: 

CgA •^'^=0 5gA -^^=0 

To prove (g^]), let s G [0, 1] be fixed; then 

Y.C{iC-x) + sp){C-x)= Yl Cisp-r]){-7j), 
CeA ri&(x+A) 

were we substituted t/ = — — x), and hence sum over —{A — x) = (x -|- A). Employing again 
(j4.9p with v{x') = —x', we obtain 

YdiC - x) + sp){C - x) = -sp, 

?ga 

and integrating with respect to s gives 

5GA •^*=0 5gA •^^=0 

Proof of Theorem 'J73 Throughout this proof, let := conv{x,^, ^ -|- and 

l|V^^i||Loc(j,^ J and ^g,^ := ||V^u||icx,(^^ j. 

We will suppress all arguments where it is possible to do so without confusion, for example, 
= S'^{u;x) and S'^ = S''{u;x). 
Defining the symbols 

Vp:=Vp{x):=Vp{Vu{x)-A,), for p G A^, j = l,2. 

we can rewrite S'^ = J2peA, ^p'^S' P- In (|4.7p we have established that ^^^\X^,pix) = 1, which 
implies 



pGA. 5eA ^ 

= E E { ['^e.p - ^p] «^ p}xU^)- (4.10) 

pGA* 5gA 

Since C has compact support, there exists a constant c > such that X^,p{x) = for all ^ G A 
with — 2;| > c\p\. Hence, we will assume throughout the rest of the proof that \C — x\ ^ \p\. 
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We Taylor expand the term [$g,p — Vp] as follows: 

^(,p-yp= 5] V^p,[-,A^^(0-V,n(x)] + i?i, (4.11) 

<;eA* 

where \E,\ < J] (|r| + |,| + 

The details of the estimates for the remainder Ei are easily established, the key observation 
being that 

\D,uiO - V,n(x)| < + le - xme^,, < (|,| + \p\)\,\ e^,„ 

which is obtained by expanding along the segments conv{x,^} and conv{^,^ + <j}. Here, and 
in the following, we skip the details required for estimating the remainders. 
Expanding D^u{S^) — V<^n(x), 

= V^_^V,u{x) + ^V^^x) + E2, (4.12) 
where \E2\ <\c^m + \p\f 6^,„ 
and combining ()4.12p with ()4.1ip yields 

- = 2^ [ • ' ^i^xVMx) + ^V2n(x)] +Ei + Es, (4.13) 

fSA, 

where < Y, '^\E.\ < ^ ^(M + Mf h., 

which we insert into (j4.10p to obtain 

- = Z E Z {^P^ [ • ' ^5-^^^ + 1^'"] ^ p}xU^) + E4, (4.14) 
pGA, geAfGA, 

where ^4 = [(^i + ^3) p]xc,p(a;). 

pSA. gGA 

Rearranging the sums, we arrive at the expression 

S'-S^= E E{^P^[-'S5eA^5.p(^)(^S-^^^ + ^^'^)]®/'}+^4. (4.15) 

peA* <?gA* 

Using ()4.7p we see that 

^Xc,pV2n = V^n, and Y,Xuy^-xy,u = -^VpV,u. (4.16) 

5eA 5eA 

Combining ()4.16p with ()4.15p . we arrive at the identity 

- = ^ E E ^P^ [ • ' ^'2/ - V, Vpy] (^p + E^. 

peA, fsA, 

Applying the symmetry Vp^ = V-p-^ (cf. ^2.3.3p . yields 

Vp,[- , V^y - V.Vpy] ® p + [ • , Vl^y - V_,V_py] ® (-/>) = 0, 

and hence we deduce that S'* — S'^ = £'4. 
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It remains to bound E4. Using its definition (j4.14p . and the bounds (|4.1ip for Ei and (|4.13p 

for E3, and the estimate Yl^£AX(,pf^ — "^^■'^CeA.x^ p(a;)7^o ^e estimate 

\Ei\ < Yl \p\Y.^\Ei\ + \E^\)x^,p 
p€A, CeA 

, ?6A 
+ ^ (|p| + kl)'"i(P,<?) max 5^,,. (4.17) 

There exists a constant c > such that, if x^,p{x) / 0) then |^ + tp — x| < c for some t G [0, 1] 
and one readily checks that this imphes 

v^^^ - X C Vp^c, ■■= {y G M"' I dist(y, conv{±2p ± q}) < 2c}. 

Hence, we obtain 

max e^.^ej.r < ||V\||ioo(3,+;,^_j||V\||ioo(j,^_^) and max (^5,? < ||V^u||2,cx>(^+;,^_^). 
xe,p(3:)7^o xe,p(a^)7^o 

Inserting these bounds into ()4.17p . and recalhng that — S'^ = i?4, we obtain the stated 
estimate. 

The statements about the sets Up^^ are easy to estabhsh. □ 

4.3. Global modeling error estimate. We present a modeling error estimate that is a nat- 
ural corollary of Theorem 14.31 For our subsequent analysis we only require the case p = 2, 
however, we give a more general statement since the same proof applies verbatim for general 
p. Earlier results this direction have been obtained in [21^ 110]. 

Lemma 4.5. Let u G W^''^ r\ K, p ^ (1, 00], and let u := ( * u; then, for all v G W^'P' , 

\{5E%u),v) - {6E%u),v)\ < (Mfp) II V^nlliP + Mi3.P) II V^ulli^,) II V^;||^,,, (4.18) 
where the constants Mg'''^^ are defined in ^2.&.4 ■ 



Proof. The case p = 00 follows immediately from Theorem l4.3l hence we assume that p G (1, 00). 
Further, we assume that v G #05 and apply a density argument to obtain the general statement. 
First, we need to show that w := G J^. To that end, we estimate 



\DpHO\ 



[ Vpu(C + tp) dt 



< |p|||Vn||Lo= < |/o|||C||li I|Vu||l°o. 



Since u G K and ||C||li = /Cdx = 1, we obtain that |DpW;(^)| < k\p\ and hence w = u\\ G Jf. 
Hence it follows that the first variations on the left-hand side are well-defined. 
From Proposition 14.11 it follows that 

\{6E^{u),v) - {5E''{u),v) \ = [ [S^{u) - S"(u)] : Vvdx 

<||S^(^i)-S^(^z)||^J|Vt-||^,,. (4.19) 

We apply the triangle inequality, 

\\S^{u) - S'{u)\\^^ < \\S%u) - S^(n)||^, + \\S'{u) - S'{n)\\^^, (4.20) 

to separately estimate the two terms on the right-hand side. 

Applying (|3.5p and Lemma lA.21 we can bound the second term by 

||S'=(Vtt) -S"(Vm)||^p < M(2)||Vn- VuIIlp < m(2)||v\||lp. (4.21) 
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To estimate the first term on the right-hand side of ()4.20p we first apply Theorem 14.31 to 
obtain 



(4.22) 



P,?eA, 



P,?,TeA, 

Next, we estimate the L^-norm of the first term on the right-hand side. To that end, we 
first recah the definition of i^p^c; from Theorem 14.31 as well as the enlarged sets defined in 

Lemma [A. 41 Let {wpi;) € ^^(A^), and let w := (Xlp^GA, w;pfmi°°'*(/>, <^))p Z^, then applying first 
Holder's inequality and then Lemma lA. 41 gives 

[ E nitHp,,)\\V^n\\L^^,^/dx<w Y w-Pm\r\p,^) I IIC*V\|1^^( )dx 



p,fGA* 



p.fSA, 



3„,I|P 
LP- 



p.fSA. 

Choosing Wp^ to balance w with ^^^^ Wp^Mp''''^ vol(fp ,-), and noting that vol(fp^) < vol( 
(this follows immediately from the definition of i^p^;) yields 



P,ceA, 

With an analogous argument we obtain 



U\\Lp- 



(4.23) 



p,?,rGA 



r 2p \ 1/p 



-7/l|2 



(4.24) 

Combining (|i:22]) . (|i:23]) and (jTM]) with ([OT]) . and noting that M^^) < Mp'^\ completes the 
proof. □ 

5. ElASTOSTATIC PROBLEMS 

In this section we present error estimates for local minimizers of the Cauchy-Born model. 
We essentially recover the result of E and Ming [W\ Thm. 2.3] for a more general class of 
interactions, and in the more challenging setting of an infinite domain and infinite interaction 
range. Moreover, due to our new consistency estimates in ^ we obtain sharper and more 
explicit estimates. 

5.1. The variational problems. 

5.1.1. Continuous external forces. For f,g & L^^^ with / • 5 € we define the inner product 

if,9)Rd := / f -gdx. 
We say that / e LJ^^ n W~^''^ if there exists a constant ||/||^^_i,2 such that 

In this case there exists a unique continuous extension of (/, • to W^''^. 
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5.1.2. The Cauchy-Born Problem. In the Cauchy-Born model, given G Lj*^^ n W we 
seek 

G avgrnin {E^{u) - {f,u)^d\ue W^'^}. (5.1) 

We understand ()5.ip as a local minimization problem with respect to the (VF^'^ H W^'°°)- 
topology. If tt'^ G -ftT is a solution to ()5.ip . then it satisfies the first-order optimality condition 

{5E%u'), v) = ir, v)^, yv G Wo. (5.2) 

We call a solution u'^ of (|5.2p stable if there exists cq > such that 

{5^E%u'')v,v) >co\\Vv\\l2 yv£Wo. (5.3) 

Prom Proposition 13.21 it follows that, if is a stable solution of ()5.2p . then it is a strict 
Pi VKi:°o)-iocal minimizer of — (f^, ■ )^d, and hence a solution of (jS.ip . 

5.1.3. External forces in the atomistic problem. Por f,gGW with f • g G we define the inner 
product 

We say that / G W^^'^ if / G #' and there exists a constant ||/||^_i,2 such that 

{f,v)A<\\f\\^^,,,\\Vvh2 yveWo. 
In this case there exists a unique continuous extension of (/, • )a to 

5.1.4. The atomistic problem. Given f'^ G we seek 

G argmin{^^(u) - {P,u)a \ u G /T^'^}. (5.4) 

We understand (|5.4p as a /oca/ minimization problem. If G is a solution to (|5.4p . then it 
satisfies the first-order optimality condition 

{5E-{u^),v) = (r, v)k yv G Wo. (5.5) 

We call a solution n*^ of ()5.5p stable if there exists cq > such that 

(52^^(n^)^;,^,) > co||V?;|||2 ^vgWo. (5.6) 

Prom Proposition 12.81 it follows that, if is a stable solution of (|5.5p then is a strict 
y^'^'^-local minimizer of i?*^ — (f '"^, ■ )a and hence a solution of (|5.4p . 

5.2. Stability of small displacements. We say that the lattice A • A is stable if 

7:= inf (5^E''{0)v,v) >0. (5.7) 

l|Vt;||^2=l 

Physically, (|5.7p states that small distortions of the Lattice A • A increase its energy. 

Por simple interactions ()5.7p can be proven analytically \26\ |9] . In practise one checks this 
stability condition by block-diagonalising S'^E'^{0) using Pouriers series l32l [lOl [15] , that is, one 
checks whether the dispersion relation satisfies uj{k) > c\k\. The condition is discussed in more 
detail in [15] and in Appendix O 

Assuming only ()5.7p we can deduce stability of "small" displacements both in the atomistic 
and Cauchy-Born models. The factor ^ in the following result is arbitrary and may be replaced 
with any number between zero and one. 

Proposition 5.1. Let A ■ A be stable, then there exists ki > such that, for k < ki, 
{5'^E''{u)v,v)>\-i\\Vv\\\2 Vz;G#^'^ Vu G ^, and 
{5'^E^{u)v,v) > ^7||V^;|||2 Vz; G W^'^, Vu G K. 
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Before we prove Proposition 15 . II we state a variant of a classical intermediate result (see, e.g., 
p. 89]; the following proof is adapted from [151 Thm. 3.1]). 



Lemma 5.2. Let 7 be defined by (|5.7|) . then 



Proof. Fix V S Wq \ {0} and set vn{x) := Nv{N ^x) for any € N, and let wn '■= vn\a', then 
we have 

{d'^E^{0)v,v) = N''^{5^E^{0)vN,VN)■ 
Takmg into account |151 Remark 1.1.1], and using the fact that v is smooth, Lemma 3.2 in |15] 
yields 

N-'^\{6^E''{0)vn,vn) - {6^E''{0)wn,wn)\ ^0 as iV ^ 00. 

We remark, that [151 Lemma 3.2] is formulated for finite-range interactions only, however, under 
the assumption that M^^^ is finite a straightforward approximation argument extends it to the 
present case. 

Using the smoothness of v it is also easy to see that N~'^^'^\\\/wn\\l^ ~^ II^'^^IIl^. Hence, we 
obtain 

^ {6^E^{0)wn,wn) n-^oo {5^E%0)v,v) 
l|V«;^||i. ^ ||V^;|12, ■ 

Taking the infimum over all v G Wq yields the stated result. □ 

Proof of Proposition \5.1[ We first consider the atomistic case. A simple variation of the proof 
of ()2.15p with i = 3 (replacing ||Vu||loo with max^gA maXpgA* \Dpu{(,)\) gives the Lipschitz 
bound 

\{{S^E''{u) -6'^E^{0))v,v)\ < cM'-^^ K\\Vv\\l2 € 

where c is a generic positive constant. Hence, choosing k < ^/{2cM^^^) yields the atomistic 
stability result. 

After employing Lemma 15.21 the proof for the continuous case is analogous. □ 



Remark 5.3. We have shown that stability of the atomistic model implies stability of the 
continuum model, using only pointwise convergence of the atomistic hessian to the continuum 
hessian (this is in fact a consequence of convergence of the energy) and scale-invariance of the 
continuum limit. Conversely, one can construct examples [3 US] where the continuum limit is 
stable (in ID, convex) while the atomistic model is not stable in the sense of (|5.7p . In this case, 
we would still expect that atomistic solutions of both the static and dynamic problem exist (in 
a suitable extended framework), however, we can no longer expect them to be "close" to the 
solutions of the Cauchy-Born equations. We give a more detailed discussion in Appendix [Ul 
from which can conclude that (|5.7p (or a similar assumption) is also necessary to obtain the 
results we seek. □ 
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5.3. Main result. We first restate the Cauchy-Born equation ()5.ip at a macroscopic scale 

X = £x, U = eu, and = e'^T, (5.8) 

where X G R"^, U, F"" : R'^ ^ R'^, and e is the atomic spacing in the X-scale. In these macro- 
scopic variables, the Cauchy-Born equation ()5.2p reads, formally 

- diYxS'iVxU') = F^ (5.9) 

By assuming that F'^ is small, more precisely, that 

II^1Ih/-i,2 + l|VF^||i2 =: 5 (5.10) 

is sufficiently small, we will be able to prove that there exists a solution to ()5.9p . Reversing 
the scaling ()5.8p . we obtain a solution ti'^(x) := e~^U^{ex) of the atomic scale equation ()5.2p . 
with external force /'^(x) := eF^{ex). We note that 

llvVlli. + llV^n^lli, =e'-'/m\Vj,U'^h. + ||V|C/^|li.) , (5.11) 

which implies that the atomistic/Cauchy-Born modelling error in the internal forces is of order 
0(e2-'^/2 (cf. Lemma To ensure that the modelling error in the external forces is of the 
same order of magnitude, we shall assume that 

\ir,v)^.-ir,v)A\<Cf6e^-''/\ (5.12) 

As a concrete example, we show in Lemma I A. 51 that, if F^ € W~^''^ H W^''^, is defined by 
dSSl), and is defined by /^(C) := / C{x - OTix) dx, then G W''^'^ and (f5l^ holds. 



(2 2) (3 2) 

Theorem 5.4. Let d < 3,k > 4, and suppose that A ■ A is stable and that Ms ' and Ms ' 
are finite (cf. § \2.3.4\ j. 

There exist constants 5o,eo > such that, for F'^ € W~^''^nW^''^ satisfying (j5.10p . /'^ defined 
by (j5.8p and G ^'-i'^ satisfying (|5.12p . and /or (5 < (5o ande < eo; ^/^ere exisi stable solutions 
u'^ and of, respectively, (|5.2p and (j5.5p . suc/i i/iai 

e<i/2||vn^- V/n^||i2 < C— , 

7 

«;/iere C = CfC{M^'^'^^/j,M^^'^^/-f). 

Remark 5.5. 1. Formally, Theorem 15 .41 states that, if the external forces are sufficiently small 
and of a "macroscopic nature" (encoded in the assumption that f''^{x) = eF'^^ex), which implies 
that V/^ V/'^ = 0(e^)), then the atomistic solution may be approximated to second-order 
accuracy by a solution of the Cauchy-Born model. 

2. The conclusion of Theorem 15.41 mav also be stated as 

\\Vu^-VIu^L^ < (Mp)||vV||L2 +Mi3'2)||v2u=||2,), 
where C" = C'{6q). In macroscopic units, with U'^{X) := eIu'^{e~'^X) the estimate reads 
live/" - VC/^||i2 < C'^e^{Mi^''^^\\V^U^\\L2+M^^''^^\\V'^U^\\li). 

On the right-hand sides of both of these estimates we may also replace with lu^, respectively 
with [/^, which effectively turns them into a priori error estimates. 

3. The factor 5 in these estimates shows that the error is O(e^) relative to the magnitude of 
the external force and hence the displacement; that is, our estimates are in fact relative error 
estimates. □ 
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The proof of this result uses a quantitative version of the inverse function theorem. The 
fohowing version is taken from [24:\ Thm. 2.1]. 

Lemma 5.6 (Inverse Function Theorem). Let si ^SS he Banach spaces, G an open subset 
of £/ , and let ^ : G ^ ^ he Frechet differentiahle. Suppose also that there exist r],a > and 
a monotone function oj : [0, +co) — )• [0, +cxd] such that 

B^{0,27]a) C G, W^my < V, < 

\\^'{U) - ,'^'{V)\\L(,^,aj) < uj{\\U - V\U) for < 27?a, 

2auj{2r]a) < 1, and auj{2rja) < 1. 

Then, there exists a unique U (z £/ such that ^{U) = and \\U\\s^ < 2r]a. 

Proof. The result follows from \2^ Thm. 2.1], upon replacing a with a~^, taking a)(t) = tuj{t) 
(admissible since co is monotone), and R = 2r]a. □ 

Proof of Theorem \5.4\ Part 1: Existence of a Cauchy-Born solution: Since A • A is stable. 
Lemma |5 . 1 1 implies that 8'^E^{Q) is positive definite, which is equivalent to the statement that 
V^Vr(O) satisfies the strong Legendre-Hadamard condition. Under this condition it is proven 
in [27] that 5^ii^'^(0) : W^'"^ — >■ W^'"^ H W~^''^ is an isomorphism. Hence, we can hope to apply 
Lemma 15.61 with 

^:=VF3,2^ ^1,2^^-1,2^ 

^{U) := 6E'{U) - F^, and G := {U £ £/ \\\U\\^ya,2 < 6}. 

Since d < 3, is embedded in VK^'°° and hence, for 6 sufficiently small, we have d K. 
Hence, (|3.5p implies that, for U,V (z G, 

\\^{U + V) - ^{U) - 6^E^U)[V, ■ ]\\yy-i.2 < M(3)||vy||^2 < M(3)||y|[^. 

A tedious but straightforward computation also shows that 

\\^{U + V)-nU)-5'E%U)[V, •]||^i,, 

= II Vdiv[VVF(VC/ + VV) - VW{VU) - V^W{VU) : VV] ||^2 

^ oidl^llj:/) for U £ G, and for sufficiently small, 

where oi{t) ^ i as t — )■ 0. (The function oi{t) depends on M^^^ and M^^^ and on the modulus 
of continuity of V^W in the set {F | |F| < k}; if A; = 5, then oi(t) = Yl%3 M'-jh^.) This shows 
that is Frechet differentiahle and 5^{U) = 5'^E''{U). 
Similarly, we can also show that 

¥^{U) - 6^{U')hi^,.j) < ooiWU - U'\U), 

where oo(t) — > as t — )• 0. (In fact, oo{t) = oi{t)/t.) 

We have in particular established that J^(0) = 6'^E^[0), which we already know to be an 
isomorphism from ^ to Moreover, by assumption we have 

\\^{0)y = \\F'\\^.,,,^^,,,<6. 

Hence, Lemma 15.61 guarantees that, for 6 sufficiently small, there exists € W^''^ H W^''^ 
satisfying (j5.9p in the strong sense, and 

IIC/'II^ < c<5/7, (5.13) 
where c is a generic constant. (The factor I/7 is due to the fact that ||5^(0)~^ ^ 1/7-) 

Let ti'^(x) := s~^U'^{ex), then the arguments given before the statement of the theorem, and 
Proposition 15. H show that is a stable solution of (|5.2p with given by ()5.8p . 
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Upon noting that ||V3;7=||^2 + WV^U^W^ < ||f/=|U ^ and that \\VU^\\l^ < \\U^\\.^ < 
5/7, we obtain the bounds 

||V=^u=||i2 + < £2-^/2^/7, and ||Vn=||Loc < 6/j. (5.14) 

Part 2: existence of an atomistic solution. Recall the definition of u'^ := * from 
Lemma SSI We apply Lemma [52] with £/ = W^'"^, = #-^'2, G ■= {w ^ #1-2 | ||Vi/;||2,2 < 
^i} for some constant 5i > that remains to be chosen. If b\ is chosen sufficiently small, then 
v^\k^ C W^"^ n J^^, hence we can define 

{^{w),v) ■= (bE''(u'' ^w),v) -{f,v)K foru;G^, v^W^'"^. 

By Theorem [221 ^ is Frechet differentiable in ^, and is Lipschitz continuous in ^, that 
is we can choose oj{t) = cM^^H in Lemma |5.6[ 

To obtain a stability estimate, we use (j5.14p and Proposition 15. II to deduce that, if 6 and 61 
are chosen sufficiently small, then 

{6'^E'''{u'')v,v) > ^7||V7;||^2 G W^'^, 

that is, ||<5^(0)-1|Il(^X) < (h)"' =■■ ^- 

To obtain a residual bound, we apply Lemmas 14.51 and IA.51 to estimate 

= {{5E^{u^),v) - {5E'^{u'^),v)] - {{r,v)A - {r,v)^,} 
< (MP)||VV||^2 +Mf 2)||vV||i. + \\Vrh.)\\Vvh.. 

that is, 

11^(0)11^ <V-=C[1 + 7-'(Mf 2) + 2))] fe2-d/2^ 
Lemma 15.61 states that, if 

cM^fl.^^-'"'^"-'"')W^<l. (5.15, 

7 V 7/7 

then there exists a locally unique solution v of c^{v) = 0. This can be guaranteed provided 
that 5e2-'='/2/7 < sq = eo{Mp'^^ /j, M^^'^^ /-f). (RecaU that M^^) < M^'^'^l) 

Let := u'^Ia- Setting n'^(^) := w'^(C) + and applying ()2.10p we obtain the estimate 
||Vn^ - < 2(77? < [1 + i(Mi2'2) + Aff'2))] l^2-d/2_ 

Applying the interpolation error estimate given in Corollary I A. 31 and ()2.10|) . we obtain 

II Vn'^ - V/n^||i2 < ||V/(u;" - 'U^)||i2 + ||Vn" - V/w;"||i2 

<||VK-n'^)||,.2 + ||VV||i2 
<[l + i(M(2.2)+Mp))]i,2-d/2_ 

This concludes the proof of the theorem. □ 

6. Convergence to solutions of the Wave Equation 
In this section we consider the dynamic problem 

{il^{t),v)A + {6E^{u''{t)),v) =0 yveWo, t>0, 
u^(0) = ng, ti^(O) = uf. 

We will prove that, if the initial condition is "macroscopic", then there exists a unique solution 
to (j6.ip . which remains close to a solution of the corresponding Cauchy-Born wave equation 
for a "macroscopic time interval". 



(6.1) 
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For simplicity, we do not consider external forces in the the dynamic problem. 

6.1. The macroscopic wave equation. Formally, the continuum limit of (|6.1|) is the Cauchy- 
Born wave equation 

u=(t) -divS"(Vu^(t)) = 0, 
subject to initial conditions. Upon rescaling 

X := ex, U := eu, T := et, (6.2) 

we formally obtain 

^U'{T) - divxS'VxU^T)) = 0, (6.3) 

which we supply with the initial condition 

[/^(O) = U^, and VtU%0) = U^. (6.4) 

To establish well-posedness of ()6.3p . ()6.4p . we apply the well-established theory. In our 
context. Theorem III in [16] reads as follows. Note that, from here on, we employ again the 
standard Sobolev spaces instead of homogeneous Sobolev spaces. 

Proposition 6.1. Let d € {1,2,3}, k > A, and suppose that A ■ A is stable. Let Uq € 
G VF3,2 ^^^f^ ||VC/^||loo < K, then there exists T'' > such that the system ([O]), ([0]) 
has a unique solution e C^{[Q,T^];W^^'^) n C^([0, T^]; 1^3,2) n C7([0, T^]; W^^'^), satisfying 
max7-g[o,T':] ||VC/'^(T)||l<» < k. 

Proof. The symbol Vt, and the conditions (al), (a2), (a3), (3.1), (3.2) in this proof refer to |16| . 
In the notation of [TU], (|6.3p reads 

a,f)=l ^ 

where aoo = L,a^p = {C{^{VU))(^^„ where C(F) = V^,W{F), h = El/3j=iSfl^' 
tio,i = fli,o = 0) for i = 1, . . . ,d. Condition (al) is trivially satisfied and condition (a2) follows 
from the fact that W £ C^. Condition (a3) is the Legendre-Hadamard condition for C(F), 
which we know from Proposition 15 . 1 1 to be satisfied for |F| < k. Hence, choosing = M"' x M'^ x 
B^dxd{0, k) we obtain (a3). 

Condition (3.1) is satisfied for s = 3. Condition (3.2) is satisfied since we assumed that 
< The requirement that s > d/2 + 1 holds since we have restricted d <2>. This 
ensures the existence of a solution with the stated regularity. 

Since G ^^([O, T^]; VF3,2^ ^ ^^([O, T^]; W^i'°°) and since ||VC/=(0)||loo < k it follows that 
||V[/^(T)||ioo < K for sufficiently short time. Thus choosing sufficiently small, we obtain 
that max^^gjo^T'c] ||Vf7^(T)||ioo < k. □ 

Upon reverting the scaling (|6.2p . we obtain the existence of a trajectory m'^(x, t) := e~^U'^ {ex , et) , 
defined for x € W^,t S [0,^^^], where := T'^/e, satisfying the following conditions: 

{u^,u)^d + {5E\u^),u) =Q VnGl^^'2, tG(0,f], (6.5) 

u^'ix^Q) = e~^U^{ex), if{x,Q>) = Ul{ex), and (6.6) 

l|V^V^tx"|lioc([o,jc]^i2) < C7„^«e'"+'"^-'^/2 for m,/ G N,l < m + Z < 4, (6.7) 

where Cm,i = \\V'^VipU''\\L^(]^Q^T'^^^L2). 
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We also imposed in Proposition 16 . 1 1 that is chosen sufficiently small to ensure that U^{T) € 
K for all T € [0, T'^]. This implies that n^(t) € K for all t G [0, t^], and hence we may conclude 
that 

{6'^E''{u^{t))v,v) > \l\\Vv\\l2 yv £f, t£ [0,tl]. 

6.2. Main result. After the preparation we can state our result on the convergence of solutions 
of (fUTT]) . 

Theorem 6.2. Suppose that d {1,2,3}, k > 4, and A ■ A is stable. Let \\VUq\\loo < k, where 

K is chosen sufficiently small so that Proposition I5.il applies. Finally, suppose that M^"''^^ is 
finite for j = 2, 3, 4, and let Cmit > be a fixed constant. 

Then there exists = £~^t^ € (OjT*^] and eo > such that, for all e < £q, and for any 
initial data Uq, nf with 

e'^/^ll v/ug - V'u"(0)||^2 + e'^'^VIul - Vn^(0)||^2 < anite^ (6.8) 

there exists G C^([0, t'^]; satisfying (|6.ip . and 

max e'^/2(||V/u^(t) - V'u=(t)||i2 + - if{t)\\L2) < Ce^, (6.9) 

where C = C{^, (M^^''))4^2, ^7^ Qnit)- 

If d <2 or k > 5, then we may choose t^ = f^. 

Proof. We give here an outline of the proof, but will establish the key technical results in the 
following subsections. All constants in this proof may depend on any property of U^. 

1. Setup: From Lemma E3] we obtain local existence for the atomistic problem: there exists 
> and li G C'^{[Q,ti];f) satisfying ([63]). Let z{t) := n'=(t)|A = (C * u%t))\A and let 

e(e,t) ■.= u^{^,t)-z{^,t), 

whenever the left-hand side is well-defined. Moreover, let w such that w{(,,t) = e{^,t). 
It follows from ()6.7p . Lemma lA.21 and Lemma \A.3\ that the initial error satisfies 

l|e(0)||,22 + i||Ve(0)||i. <Coe^-^ (6.10) 

where Cq depends on Uq, Uf, Cinit but is independent of e. We fix a constant > Cq, which 
will be specified later on. Since u'^ G C^([0, ti], ^^) (and hence e G C^([0, ti], ^^)) upon choosing 
ti sufficiently small, we obtain 

l|e(*)ll'2 + il|Ve(t)||22 < Ce^-"^ Vt G [0,ti]. (6.11) 

Thus the task is to prove that we may choose ti > e^^ . 

From the definition of z and the assumption that || Vn^(t) |[l=o < k! < k for all t G [O,^'^], we 
deduce that \Dpz{t,^)\ < k' for all t,^,p (cf. the proof of Lemma 113]) . Choosing eq sufficiently 
small, then for e < Eq we obtain from ()6.1ip that 

z{t) + ee{t) e VtG[0,ti], 6Ig[0, 1]. (6.12) 

Since we assumed that k, was chosen sufficiently small. Proposition 15.11 implies that 

{6^E''{z{t)+ee{t))v,v) >^Vv\\l2 G ^^ t £ [0,ti], 0g[O,1]. (6.13) 

2. Error equation: Testing ()6.ip with e and ()6.5p with w, yields 

{e,e)A + {6E^{u^)-5E^{z),e) = {{il',w)^, - (z,e)A} + {{6E%u'),w) - {6E^{z),e)}. (6.14) 
Since w = e, the first group on the right-hand side can be rewritten as (cf. (]A.3P ) 

[iu'','w)Kd - iz,e)A^ = (u'' - z,w)^d =: {a,w)^d. 
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where z and w are identified with their first-order interpolants. Similarly, using Proposition 14. 1|, 
we can write 

- (<5i?'^(z),e)} = (S=K) -S'^(z),Vu;V =: (^,Vu;V, 
Moreover, according to ()6.12p we can expand 

((5E^(n^) -(5E''(z),e> =: / (^^^^(z + 0e)e, e>d0 =: (Fe, e), 
to rewrite (j6.14p as 

-^{||e||i + (i7e,e)} = (i^e, e) + 2(a, e)iR. + 2(/3, Vu;)^^. (6.15) 
We define E'^{t) := ||e||^2 + {He,e), and note that ()6.13p immediately implies that 

\\e\\% + i\\ye\\l2 <E'^{t) <CH{\\e\\% + ^\\Ve\\l2) VtG[0,ti], (6.16) 
where Ch depends only on M^'^\ We integrate (|6.15p over (0,s), s < ti, to obtain 

E^{s) = E^{0) + \^{He,e) +2{a,w)^d + 2{/3,w)^djdt, (6.17) 

and proceed to estimate the three terms on the right-hand side separately. 

3. Consistency estimates and nonlinearity: The second and third term on the right-hand 
side of (|6.17p measure the consistency between the atomistic and Cauchy-Born model. Stan- 
dard interpolation error results (cf. ^6.3l for the details) yield the estimate ||a||L2 < ||V^ii'^||i2 < 
g.3-d/2^ and hence, applying Cauchy's inequality and (j2.5p . we obtain 

^ 2{a,w)^ddt <Ca(^tie^-'^ + e E'^dt^, (6.18) 



where Ca depends only on the trajectory f/'^ through the bounds (|6.7p . but it is independent 
of £. 

The third term on the right-hand side of ()6.17p requires an estimate that it similar to Lemma 
We establish the required variant in 



w\\l2, (6.19) 

where C = CiU", (M^^''^^)^^2)- Integrating over (0, s), s < ti, and applying Cauchy's inequality, 
we again obtain 



Jo 



2{I3,\/w)l2 dt < C^(^ti£^-'^ + £ E'^dt^ (6.20) 



where depends on (Mi'''^^)^^2 ^^d on W^, but not on e (provided £ < £o, which we chose 
above so that = z + e ^ J^) . 

Finally, the first term on the right-hand of ()6.17p side is not a consistency term, but must be 
otherwise controlled. Our argument is a refinement of a method due to Makridakis [22] for the 
numerical analysis of nonlinear wave equations. In i ^6.4l we show that an integrating by parts 
argument leads to 

^ {He,e)dt<C,,i(^E^{s) + E^{0) + (^£E^ + £E^ + E"^^ dt^ , (6.21) 
where C^i = C^iiM^^^ M^^^U"). 
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4. Gronwall lemma: Combining ([6l8]) . ([OO]) . ([OT]) with I^J7\) yields 

E\s) < E\0) + {Ca + Cp) (tit'-^ + ^ 

We employ (l630]) . (I6TT]) . (1636]) and Ti := tie, to obtain 

E\s) < [ChCo + (C« + Cp)Ti] e^-'' + e(C^ + C^j + Cni) /" i?^ dt 

JO 

+ Cnl (c^/'Co'/'66-3d/2 ^ ^7^/2^3/2^6-3^/2 ^ C^C^e^-'^i^^ dt) . 

Since £6-2ci/2 ^ . g2 (i/2 gince d < 3, choosing eo sufficiently small ensures that 

C|/'C7o'/'e6-3d/2 ^^3/2^3/2^6-3d/2 < C'^^Coe^"'^. 

(However, if d = 3 then e^"*^ = e, hence the integral term cannot be made arbitrarily small 
for s = 0(l/e).) Upon defining Ci = Ci{Ti) = 2ChCo + {C^ + Cp)Ti and C2 = C2(a,eo) = 
Co. + Cp + Cni(l + ChC.eI-''), yields 

E'^{s) < Cie^-'^ + C2 [ E^ dt, for < s < ti. (6.22) 
Jo 

Applying Gronwall's inequality, we obtain 

max E^{t) < exp(C72ri)e^-'^. (6.23) 

We observe that Ci exp(C2ri) 2ChCo as Ti 0. 

We now choose C^,eo,T''^ in the following order: 1. := AChCq (or any constant larger 
than 2ChCo); 2. eo < 1 and sufficiently small as required in the estimates up to this point 
(dependent in particular on C*); and 3. > sufficiently small so that 

Ci(r^) exp(C2(C7,, l)r^) =: CsiT^) < C, and < 

This is possible due to the fact that C3(r^) as 0. We set := T^/e. 

Using these definitions, we can estimate (|6.23p above by 

max E\t) < C3{Ti)e^~'^ < C,e^~'^. (6.24) 

5. Continuation argument: Suppose that ti < t^ is chosen maximally, such that the 
trajectory exists in [0, ti] and (|6.1ip holds with the choice of C^, we made above. If ti < t^, 
then we end up with the stronger error estimate (|6.24p . We apply Lemma 16.31 again to extend 
u'^ to an interval [ti,t2] where ti < t2 < t^. Since ()6.24p holds in [0, ti], and since Cs{Ti) < 
when Ti < T^, choosing t2 sufficiently close to ti implies that (|6.1ip holds in [0,t2]- This 
contradicts the maximality of ti and hence we must have ti = t^. Thus, we conclude that 
(IHTB holds in [0,t^], where = e'^T'^ and > 0. 

6. Maximal time interval for d < 3: If c? < 3 then one varies the argument following ()6.23p . 
exploiting the fact that e^"'^ ^ as eo ^ 0, as follows. 1. Choose C* := Ci(T'=) exp(C2(l, l)r'=). 
2. Choose eo sufficiently small as required by the previous estimates, and in addition sufficiently 
small so that C,,eo < 1- 3. Choose = T^. The rest of the argument is analogous. 

7. Maximal time interval for k > 5: Now suppose that d = 3 and k > 5, then we need to 
vary the proof starting from (|6.2ip . The assumption that k > 5 allows us to integrate the term 
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{He,e) by parts twice, to prove in ^6.41 that 

{He, e) dt < Cni(^E=^(s) + ^^(o) + E^(s) + E^(0) 



where Cni depends only on (M^^^)^^^. Exploiting the higher powers of E and arguing similarly 
as in 4, upon choosing £q sufficiently small, we obtain 

[He, e) dt < ChCqe^''^ + C,,xe [ E^ dt, 

Jo 

and eventually 

max E^{t) < Ciexp(C72Ti)e^"'^, 

where Ci = 2ChCq + {Ca + Ci3)Ti and C2 = Ca + Cp + Ca\. Since Ci and C2 are independent 
of C* we can now argue as in step 6 to deduce that we may choose = T'^_ □ 

Lemma 6.3. Let uq,ui S £'^, ||Viio||L°o < k, then there exists ti > and a unique trajectory 
u € C'^{[0,ti];f) satisfying (f6TT]l with n(0) = Uo,u{0) = Ui. 

Proof. The estimate (j2.15p with j = 2 implies the existence of 5 > such that, for all vi,V2 € 
Be2{uo,d), 

lidE^'ivi) - 5E^{v2),w)\ < M^^^llVt;! - Vv2\\l2\\Vw\\l2 < M^'^^\\vi - 11^2 ||'u;||^2. 
that is, dE^ : ^ ^ is Lipschitz. Hence the result follows from Picard's theorem. □ 

6.3. Consistency estimates. In this section, we prove the two consistency estimates ()6.18p 
and (|6.20p . We begin by considering ()6.18p . which follows from standard interpolation error 
estimates. 



Lemma 6.4. Let v € W'^''^, and let C,*v he the first-order interpolant of Q * f |a, then 
Proof. We first apply a triangle inequality. 



L2- 



The first term on the right-hand side is estimated using Lemma \A.2\ 

\\'"-C*v\\l2 ^ l|V^t^|lL2. 

The second term on the right-hand side is a standard Ql-interpolation error, which can be 
estimated by [5] 



\\C*v-C*v\\^^ < \\V\C*v)\\l2 = lie* V'7;|Ii2 < \\V'v\\l2. □ 

Proof of ()6.18p . Applying Lemma 16.41 the regularity bound (|6.7p . and the norm-equivalence 
II^IIl2 ~ 11^11^2 < E (see S 12. ip . we can estimate 

f S f s p s 

/ {a,w)L2dt< / ||VV||i2||u;||i2dt < / e3-d/2^dt. 
Jo Jo Jo 

Applying Cauchy's inequality to e^~'^^'^E = e(^~<^)/2 . (e^/'^E), we obtain 

^ (a, w)l2 dt< (e^""^ + e^^) dt < tie^~'^ + e E^ dt □ 
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We now turn towards (j6.20p . which is the crucial ingredient in the dynamic analysis. Our 
strategy will be to integrate by parts, and estimate V/3, which is the error in the divergence of 
the stresses. To that end, we first need to prove that x ^ S'^{u]x) is differentiable. For S'^, the 
analogous result follows from Lemma |3. 11 

Lemma 6.5. Let u € , then x ^ S^{u]x) e with 

divS^(n;x) = (^.26) 

pGA, geA 

where V pX^,p{^) = ^X^A^) ' P- 

Proof. Since xg.p ^ and (x5,p)ceA is a partition of unity, it is straightforward to show that 
|S'^(m; x)\ < M(i) for ah x € M*^. 

The rest of the proof is devoted to the assertion that S'^ is Lipschitz continuous. First, 
we compute the Lipschitz constant of Xi,,p- Let L be the global Lipschitz constant of let 
Sx ■■= supp(C( • - x)), and let J^^y := {s G (0, |p|) | ^ + G 5^ n Sy}; then 

{x)-xUy)\<^\ [ |C(e + S|^-x)-C(^ + s^-y)|d5 

. 1 /■ ,1 I , . 2Ldiam(^ 

Next, we note that ^^(z\{xs,,p{x) ~ X^,piy)) = 0, and hence we can rewrite 
S'^(n; x) - S^{u; y) = 2^[^?,p W ® fe,p(^) " 

pGA* gGA 
pGA, gGA 

Expanding using that fact that 

#{e e A I 7^ and x?,p(y) 7^ 0} < [pj, 

and the Lipschitz estimate for Xc,p; we obtain 

\S^{u;x) - S^{u;y)\ < M^^^\\Vu\\l^\x - y\. 

The formula for the divergence (j6.26p is easy to establish. □ 

We see from the formula for divS'^ that we will require symmetries of the weighting functions 
VpX?,P) which we establish next. 

Lemma 6.6. Let x G M"^ and p G A*, then 

Y'^pXU^)=0, (6.27) 

Y ^^pX^A^) - 2;) = p, and (6.28) 

eeA 

Yl ^pX^A^) {^-x)(^{^-x) = -p(^p. (6.29) 
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Proof. First, we compute VpX£^,p{x): 

VpXC,pW= - / VC{^ + tp-x)pdt 
Jo 



C i-Ai +tp-x)dt = cii -x)-ai+p-x) 

Jo 



.30) 



Using this identity and the fact that the nodal interpolant reproduces affine functions, we 
immediately obtain (j6.27p and (j6.28p 
To prove (j6.29p . consider 

= E - ^) - c(e + P - x)] (e - x) ® (e - x) 

= c(e - 2:) [(c - 2;) ® (e - x) - - p - x) ® (e - p - x)] 

= Y CiC - x)[p ® - x) + - x) ^ p - p ® p]. 



Using again the fact that the nodal interpolant reproduces affine functions, we obtain 

VpX5,p(x) - x) ® - x) = p (gi {x - x) + {x — x) <Si p - p ^ p = -p® p. □ 



The key step towards proving (|6.20p is the next result, which is analogous to Lemma 



Lemma 6.7. Let u € W,^'°° D K; then 



loc 

4 

(00), 



|divS-(u;x)-divS=(^;x)| < J] ni^P ip)T{\\'^ML^(.x+u,,,^) 



PGA2 



(x+l/p,,p2)ll V "|1L°°(X+I^p,,p3) 



where is defined in { 2.3.4\ 



Proof. Recall from the proof of Lemma 14.31 the definition of the sets I'g^^ and fp^,; and the 
quantities e^^^ and 5^^^. In addition, we define ^^^^ := || V^u||j;^oo(-^,^ j. 

Recall from (|6.26p the algebraic expression for divS'^. Since S'^ = X^pg^^ ® /?, where 
Vp := Vp{Vu{x) ■ A^), we also obtain 

divS^(x)= ^ yp,[. ,VpV,7x(x)]. (6.31) 

p,?ga. 
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To proceed, we expand ^^^p in (|6.26p to third order: 

feA, 



<;,r6A, 

=: Vp + ri(e, p) + p) + E2{^, p) (6.32) 
where |£;2(^,p)|< ^ "^^tl""'^ HC^' + '^-D "S-^-' 

P2,P3,P4eA, j=2 

where the estimate for £"2 is fully analogous to the estimate of Ei in (j4.1ip (taking into ac- 
count the finite support of X5,p)- Summing over ^ € A,/> G A*, employing the fact that 
^$eA W pXi,p\ ^ 2 (this is an immediate consequence of ()6.30p ). and otherwise arguing as in 
the proof of Theorem 14. 3^ we obtain 

peA*,eeA peAf iii=i l/^^UeA i=2 

^ E 11 ^1/ Ei^pi>^5,Pii 

< )(p) n l|V^^|U^(.+.,^„j. (6.33) 

Moreover, using (|6.27p . it is straightforward to treat the first term in (|6.32p : 

E Y.^pXU^)Vp = Q. (6.34) 

peA, 5gA 

We will now independently estimate the two remaining groups involving Ti and Ei. 
Estimating Ti: We expand 

D,u{i) - V,u{x) = (Vg.^V, + iV^ + |V|_,.V, + + \Vf)u{x) + E^i^c;), 

where mC,<;)\ < {\p\ + \,\f\,\p^,,. 

Summing over S^, p, applying (j6.27p - (j6.29p . and applying the symmetry Vp^ = V^p, yields 

EE^pX«.P^i(^'/')= E Vp,[ - ,iVpV,-'^VlV, + ^VpV^M^)]+E', 
psAt ^eA p,f eA, 

= E ^p.[- ,VpV,u]+E^ = divS^ + E^, (6.35) 

P.^eA. 

where {E'^l < E ^^d°\p) \\^^u\\l^{x+u,,,p^), 

and where the estimate of \E'^\ is analogous to the estimate of \E'2\ in ()6.33p . 

Estimating Ei: To estimate Ei we use the simpler expansion ()4.12p . which can be written 

as 

Dc;u{i) - V^u = Vi:_^V^u + i V^n + 5'^ 
where <{\p\ + <i\?m,,. 



2, (6-36) 
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We apply again (j6.27p - (j6.29p . followed by the symmetry Vp^r = —V-p-i-T, to obtain 
peA, 5eA p,f,reA* ^eA 

+ ^Vp^A ■ , V^n, VpVru] + • , V^n, V^n]) + E[ 

= E[, (6.37) 

where 

p,?,tgA* 5eA 

p,f ,rgA, ^gA 

P,?,tGA. 5ga 
=: £^1^1 + E[^2 + -E'i_3. 

The first and second terms can be treated in a straightforward manner, employing ()6.36|) : 

Ki\ ^ E 4"^(p)|V'n(x)j||V3^x||i^(,+,^^,^^), and 

However, the term E[ 3 must be treated more carefully. Note that simply using (|6.36p to 
estimate 

would impose a decay on the interaction potential that would rule out the Lennard-Jones 
potential. Instead, we use ()6.36p only for the terms 5'^ ^ and the weaker estimate 

14^,1 = \DrU{0 - VrU - V^^^VrU - iV^nj < (|p| + \T\)\T\e^^r, 

which leads to 

ip/ I < \- "^(P)lpl^ |y73, n ||V72„|| 

peA3 

Note that the coefficients in this last estimate are precisely m^°°^(p), which are controlled also 
for the Lennard-Jones case. Combining the estimates for E'^^,E'^2 ^'^^ -^1 3' arrive at 

1^11 ^ E 4"^(P)l|V'^IU-(x+.p„p,)l|V't^llL-(x+.,,„3). (6.38) 

Note that we have seemingly ignored the term |, however, due to the symmetry 'iTi{pi, p2, ps) = 
m{pi, p3, P2), this term is in fact included. 

Combining the estimates: Combining (j6.33p . ()6.35p and ()6.38p we obtain the desired upper 
bound on IdivS*^ — divS'^l. □ 
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Repeating the arguments of the proof of Lemma [4 . 5 1 almost verbatim, we obtain the following 
global consistency error estimate. 

Lemma 6.8. Let u G W^'^ , u := C * u\\ and v S £p' ; then 



[5E\u),v) - {5E''{u),v)\ < (Mf ^)||V^n||iP + Mf ^)||V\||i3p/2||V2n| 



L3p 



+ Mf^)||V\|li3.)||.||^,.. 

We are now in a position to prove ()6.20p . 
Proof of ([OO]) . Recall that 

{P,Vw) = {6E^{u^),w) - {6E-{z),e), 
where w = e and z = {(^ * u'^)\a- Hence, Lemma [6^ and ()6.7p yield 

\{P,Vw)\ < Ce^-'^/^\\w\\L2, 

where C depends only on Mj^''^\ j = 2,3,4, and is otherwise generic. This proves (j6.19p . from 
which (j6.20p follows immediately after application of Cauchy's inequality. □ 

6.4. Estimating the nonlinearity. In this section, we prove two estimates on the term {He, e) 
occurring in the error equation ()6.15p . 

Proof of (fOTjl . Recall (f632]) . which states that z{t) + ee{t) € JT for all t € [0, s] and 6 € [0, 1]. 

We write out {He, e) explicitly and exchange the order of integration (this is justified since 
M^^) is finite, which implies that the integrand is uniformly bounded), 

/ {He,e)dt= / / d^E^{z + ee)[z + ee,e,e]de dt 
Jo Jo Jo 

pi ps pi ps 

= 6^E''{z + ee)[i,e,e] dtde+ 5^E^{z + ee)[e,e,e]dtde. 

Jo Jo Jo Jo 



The first group on the right-hand side is easily estimated, using (|2.15p and (|6.7p . 

6^E''iz + ee)[z,e,e] < Af^^^ || Vi||Loo ||Ve||^2 < sM^^^\\VU''\\l^E^. 

To estimate the second group, we focus on the inner integral only. Since 5^E''^( ■ )[vi,V2,V3] 
is invariant under permutation of the arguments vi, V2,V3 it follows that 

-^l6^E^{z + ee)[e,e,e] = 6^E^{z + ee)[e,e,e] + y^E'^iz + ee)[z + 9e,e,e,e], 

which gives 

ps ps 

/ 5^E^{z + ee)[e,e,e]dt = ^[6^E^{z + ee)[e,e,e]Y -I / 6^E^{z + 0e)[z + ee,e,e,e]dt. 
Jo Jo 

(6.39) 

Applying (j2.15p and the inverse estimate Lemma 12.31 with j = 0,m = 1 and p = oo, g = 2 to 
bound ||Ve||Loo < ||e||i2, we obtain 



'el\h 



r d^E'^iz + ee)[e, e, e] dt < M^^) (|| Ve(0)||i2 + ||Ve(s) + M^^) /" \\Vz\\lo. || Ve 
Jo Jo 

+ m(^) [ ||Ve||Loo||Ve||i2dt 
Jo 

<M(3)(E(0)3 + ^(s)3) +m(^)^ (^eE^ + E^^dt. □ 
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Proof of (j6.25p . To prove ()6.25p we may now use the fact that k > 5, that is, i?^ is five times 
differentiable in a neighborhood of z{t) + 9e{t), 9 G [0, l],t € [0, ti]. Thus, starting from (j6.39p 
we can repeat the integration by parts argument to obtain (|6.25p . □ 

7. Conclusion 

In the small deformation or short time regimes we have developed an essentially complete 
approximation theory of the Cauchy~Born model for Bravais lattices for both static and dy- 
namic problems. The main open questions are 1. the extension to large deformations and/or 
the charactisation of maximal time intervals for which the error estimates hold; and 2. the 
extension to multi-lattices. 
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Appendix A. Auxiliary Approximation and Interpolation Results 

In this appendix we collect several auxiliary results related to the smoothing properties of 
the convolution operator f ^ Q * f , and the accuracy of the approximation C,* f to f . 

Lemma A.l. Let C, he the basis function defined in H2.1\. and f{x) = a + b ■ x where 
a G M, 6 € M'^ are constants, then Q* f = f ■ 

Proof. Since / C,{x)dx = 1, it follows that C, * a = a. Further, 

/ Cix-y)ydy=x- C{x - y){x - y) dy = X - ((-y) ■ {-y) dy = x, 
where we used that fact that ^ is an even function. □ 

Lemma A. 2. Let f G W'^''^, p G [1, oo], then 

IIC*/-/I|lp < IIVVIIlp, for j = 0,1,2. 



Proof. This result follows immediately from Lemma I A. II and Poincare inequalities. □ 



Corollary A. 3. Let u G VK^'P with p > d/3 and let u := C, *u, then 

\Vu-VIu\\lv < \\V^u\\lp. 



Proof. Under the assumption that p > d/3 (which guarantees the embedding W^'^ C C) we 
know from [25] that 

\\Vu - V/uIIlp < W^uWlp < W^Alv. 
Applying also Lemma lA. 21 with / = Vu, we obtain 

\\Vu - VLu\\lp < llVn - V-uIIlp + llVu - VIu\\lp < ||V^n||LP. □ 
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Lemma A. 4. Let u cM. be measurable, —v = v, and f € LP(M. y , then 

llC*/ll^,^,.^,.dx<vol(z.') 



where v' := []^^i,swvv{C{x - ■))■ 

Proof. Let C,'{x,z) := may^y^x+v Qiv — z). Since C ^ ™d J^dQ^x = 1, Jensen's inequality 
yields 



ay-z)f{z)dz 



< max / C{y-z)\fiz)\Pdz 



< maxC{y-z)\f{z)\Pdz= ('{x, z)\f{zW dz. 
Integrating with respect to x, we obtain 

lC*/llio.(,+,)dx< / / C'{x,z)\fiz)\Pdzdx 

" / c'{x,z)dxdz. 



From its definition it is clear that < < 1. Moreover, if C'{.x, z) ^ 0, then C,{y — z) 7^ for 
some y ^ X + V, that is, 

y — X € u and z — y £ suppC- 
Since both v and suppC are symmetric about the origin, this implies that x — z G \Jy^,y supp(y — 
• ) =: I/'. Thus, we obtain J^d C{x,z) dx < vol(i^'), which concludes the proof. □ 

Lemma A. 5. Let € Lj^^ n W~^''^ and let e W be defined by 

r(c):= / a(-x)r{x)dx, (a.i) 

then e Moreover, ifV € L^, then 

|(r,^)A - {f\v)^d\ < iivriii2||v^k2 € ^0. (A.2) 

Proof. Let v gWq, then 

{r,v)^ = Y,ni)-v{i)= I nx)Y.a^-i)v{i)dx= {r,v)^d, 

5GA ceA 

which implies that G #-^'2 with Wf'Wy^-ia < ||/11vi/-i,2. 
To prove (|A.2p . we first note that 

(r,^)A = J^ne) / c(x-eMx)dx 

«eA -^^^ 

= [ v{x).^ax-Ofmdx = {f^v)^d. (A.3) 
jRd CeA 

This allows us to write 

{r,v)A-{r,v)^. = Yl I c{(-x){f^{o-f'{x))v{x)dx. 
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From the definition of (jA.ip it follows that 

/ c(e-^)(r(e)-r(x)) = o, 

and hence we obtain, for some arbitrary constants G M, 

\{r,v)A-{r,vhA 



E / c(e - x) (r (c) - rix)) (n(x) - c^) d. 
^ E c(e - X) - r(x)p dx) ||. - c^||^,(Q^). 



Choosing = (v)qj and applying Poincare's Inequality, we can estimate 

ll^-cdL^cQ,) - (^)1^HIl2(q,)- 

Moreover, estimating C ^ li and using the fact that /^(^) is the orthogonal projection of 
with respect to the kernel C('^ ~ " )i we obtain 

<[^j c(e-x)|(r)Q,-r(x)pdxy^' 

Combining the foregoing estimates and estimating the overlaps we arrive at 

\{r,v)K-{r,v)^,\<c\\vnLAWv\\L^. 

with C = < 1 for (i G {1, 2, 3}. This establishes (ICT . □ 

Appendix B. Examples of admissible potentials 

We discuss the most common interatomic potentials and show that they can be accommo- 
dated within our framework. We remark from the outset that our smoothness requirement 
(at least four times continuously differentiable for the error analysis) is reasonable for physi- 
cal interaction potentials, but is not satisfied by typical potentials constructed for molecular 
dynamics simulations, which employ cut-off functions that are often only once differentiable. 

B.l. Lennard-Jones type potentials. For a pure pair interaction model, we define 

yi9):=lYl H\9p\)-H\^P\)]- (B.l) 
peA, 

With this definition, V clearly satisfies the symmetry (j2.1ip . 

The prototypical example is of course the Lennard-Jones potential |17j . 

In this case, one readily sees that ip^^\r) < r~^~^ , for r > 1 and j G N. More generally, suppose 
that V is of the form (jB.ip with 

< ^-"-i^ forr>l, j G N, (B.2) 

then one may readily deduce that 

l/ol"", p= (p, ...,/)) G Ai, 



m{p) 



0, otherwise, 
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and consequently, 

M^^^ < \P\~" and + M^^''^ < \p\~"+''/'^. 

Thus, M^-^^ is finite (and hence well-defined and k times differentiable) if and only if a > d; 
and M^^''^\ Mi^'^^ are finite (and hence our error analysis applies) if and only if a > d + 5/2. 

In particular, it follows that the Lennard- Jones potential is included in our analysis. Another 
commonly employed potential is the Morse potential [23], which decays exponentially and is 
hence trivially included our analysis. The Coulomb potential, ip{r) = r~^, is excluded. 

B.2. Embedded atom method. In the embedded atom method [7] one postulates site ener- 
gies of the form 

V{9)= E<^(I^'pI) + ^(Sp6aXI5pI))' 
peA, 

where ip is a. Lennard- Jones or Morse type pair potential, i/'dgpl) is a model of the electron 
density at generated by a nucleus at distance j^fpl, and G is the energy to embed a nucleus 
into a see of electrons. Again, it is clear from the functional form of V, that it satisfies the 
symmetry ()2.1ip . 

The computation of the partial derivatives is now more involved. Suppose, for simplicity, 
that ip = 0, define ip := EpeA. ^(bpl), and "^{g) := V'dfl), then 

Vp{g)=G'{i^)V^{gp), 

Vpdg) = G"{i;)V^{gp) V^{g,) + G'{i^)V^^{gp)6p,„ 
Vp,r{9) = G"'{^)V^{gp) ® V^{g,) Vf (^0 + G' {^)V^^{gp)6pJp,r 
+ G"{^) (V^^igp) ® V^{g,)Sp,r + V^igp) ® V^^{9,)S,,r + ^^"^{Op) ^ Vf (5r)<5p,,) , 

and so forth. G is typically chosen smooth and 93, ip decay exponentially. In that case, one 
immediately sees that all constants M^^), Mi^''^\ M^^'^^ are bounded. More generally, let <^ = 
and suppose that 

< (B.4) 
In this case, also \V^"^{g)\ < \g\~^~'' , and we conclude that 

m{p,^) <\p\~'^\^\~'^ + \p\~'^6p„ p,? € A,, 
Hp, r) < \p\-^\<^\-^\t\-^ + \p\-^SpJp,r 

+ (ipr^kr^v? + \p\~^\r\-^^p,r + kr^kr^,.), p,?,t g a„ 

and so forth. Due to the product structure, one can readily see that M^^^ is finite provided 
that /3 > d. 

(7 2) (j 2) 

However, to ensure that , Ms are finite, we now require more stringent requirements. 

For example, considering only the first group in m{p, <f, r) and indicating the missing terms by 
". . .", and using |/o x <f| < \p\\<;\, we can estimate 

Mi^'^)< y: (H + ki + M)^ipr^kr^irr^ + ... 

p,f,r6A, 

p,?,TeA« 
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which is finite provided that /3 > d + 3. The remaining terms can be treated analogously. For 
the dynamic case, the extra factor |/0i|~^ does not help except in the case of pair interactions, 

and we require /3 > d + 4 to ensure that the constants M^'^^ are finite. 

In summary, if V is of the form (jB.3p with the pair interaction (p satisfying ()B.2p and 
the electron density function ip satisfying ()B.4p . then we require q,/3 > d to ensure that the 

constants M^^\ 1 < j < k are finite; we require a > d + 5/2, /3 > d + 3 to ensure that Ms''''^\ 
j = 2, 3, are finite; and we require a > d + 5/2, /3 > d + 4 to ensure that M^"''^^ 2 < j < 4 are 
finite. 

B.3. Bond-angle potentials. Lennard-Jones type pair interactions and embedded atom po- 
tentials are the prototypical long-ranged potentials with infinite interaction range. Most other 
potentials used in molecular simulations act only on a finite interaction neighborhood. For 
example, bond-angle potentials (3-body or 4-body) act only on angles between nearest neigh- 
bors. We only need to check whether they can be written in a way that preserves the inversion 
symmetry (j2.1ip . 3-body bond-angle energies are typically written in the form 

where r^^ = rj + u{r]) — ^ — 'u(^), 9^^^ is the angle between the bond directions r^^,r^^, tp 
is a cut-off function to ensure that the potential acts only on nearest-neighbors, and tjj is an 
angle potential that drives towards preferred bond-angles. This term is symmetric about the 
center-atom, which suggests to write 

V{g)= v{\rp\M\r,\mep,), 
P,?eA. 

where rp := p + gp, and 6p^ is the angle between rp, r^. This sum is fully permutation invariant, 
and hence the inversion symmetry (|2.1ip holds. 

4-body (or, dihedral angle; or, torsion) potentials can be treated similarly. There are now 
two center atoms in bonds of this type, and hence one "splits" the bond between the two 
corresponding site energies (similarly as in the pair potential case). By summing over all 
quadruples involved with the given site, the resulting site potential will again be permutation 
invariant. 

B.4. Generic multi-body potentials. A more recent development are potentials without 
physical interpretation, but simply postulating a general functional form for V, and fitting a 
large number of parameters to energy and forces obtained from electronic structure calculations; 
see, e.g., [2]. Such general potentials are normally constructed to satisfy the permutation 
invariance, and hence the inversion symmetry (|2.1ip . and are therefore still included in our 
analysis. 

Appendix C. Lattice Stability versus Ellipticity 

We show that the lattice stability assumption (|5.7p is not only sufficient but also necessary 
to obtain Theorems 15.41 and 16.21 This can already be seen for ID second-neighbour harmonic 
pair interactions: 
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where = — In this case, the atomistic and Cauchy-Born energies are more conve- 

niently written in the form 

E-{u) = Y,{W^\^ + fK + 4+i\^) and E'' (u) = ^-^^^^ f \u' \^ dx . 

?gA 

We consider two choices for the coefficients 01,02: 

(1) _ _i and (2) _ 1 

"■2 ~ 4 ' "2 ~ 2 ' 

then in both of these cases we have 

a[^^ + 44^') = 1. 

Thus, the continuum energy is positive definite and hence the continuum wave equation is 
weh-posed. 

In the atomistic case, we can use the parahelogram formula to rewrite 

E^{u) = Y,{^^'^W'^\^ -f\u'l\^), (C.l) 

where u'^ = n^+i — 2u^ + u^-i- Hence, in the case = a^-^^ we have that E^{u) > E'^{u), so 
that (|5.7p is satisfied and the dynamic atomistic and continuum solutions will remain close for 
a macroscopic time interval (cf. Theorem I6.2p . 

By contrast, in the case = , where 02 > we can see from (|C.ip that oscillations are 
energetically advantageous. Indeed we note that, formally, defining (p'{(,) = (—1)^ gives infinite 
negative energy. 

Formally (since (p ^ W^''^) one easily checks that Hip = —ip, where H := 5'^E'^{0). A straight- 
forward approximation argument shows that —1 belongs to the spectrum of H. 

For the static case, this means that even if atomistic solutions exist, they are not local 
minimizers. 

For the dynamic case, it means that there exist exponentially growing solutions. Using the 
characterization of the spectrum in terms of approximate eigenfunctions, there exists for each 
S > a function tps G with HV'^Ufa = 1, such that \\Hips + ipsW ^ ^ (see Section VIII. 3 in 
[30]). Suppose now that we solve the Cauchy-Born equation with m(0) = ■it(O) = and the 
atomistic equation with ti(0) = and u{0) = e^il^^. The function v{t) := sinh(t)e^?/^5 then solves 
the atomistic evolution equation to order 0{5). By estimating the difference u{t) — v{t) it is 
straightforward to prove that 

\\u(t)\\p > e^ie* for t < 3|loge|, 
and in particular, ||ti||£2 becomes of order one for t \ loge|. 
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